Noname manuscript No. 

(will be inserted by the editor) 



Smart expansion and fast calibration for jump diffusion 

E. Benhamou E. Gobet M. Miri 



First version: May 11, 2007/ This version: September 17, 2008 



Abstract Using Malliavin calculus techniques, we derive an analytical formula for 
the price of European options, for any model including local volatility and Poisson 
jump process. We show that the accuracy of the formula depends on the smoothness 
of the payoff function. Our approach relies on an asymptotic expansion related to 
small diffusion and small jump frequency/size. Our formula has excellent accuracy 
(the error on implied Black-Scholes volatilities for call option is smaller than 2 bp for 
various strikes and maturities). Additionally, model calibration becomes very rapid. 

Keywords asymptotic expansion • Malliavin calculus • volatility skew and smile • 
small diffusion process • small jump frequency /size 

JEL Classification: G13 

Mathematics Subject Classification (2000): 60J75 60HXX 



1 Introduction 

The standard Black-Scholes formula (1973) was derived under the assumption of 
lognormal diffusion with constant volatility to price calls and puts. However, this hy- 
pothesis is unrealistic under real market conditions because we need to use different 
volatilities to equate different option strikes K and maturities T . Market data shows 
that the shape of the implied volatilities takes the form of a smile or a skew. 

In order to fit the smile or the skew, Dupire (in | |Dup94) ) and Rubinstein (in 
IIRub94l ) use a local volatiUty CT/„e(f,/) depending on time t and state / to fit the 
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market. This hypothesis is interesting for hedging because it maintains the complete- 
ness of the market. However, in a few cases IIACCLOll . one has closed formulas. In 
the case of homogeneous volatility, singular perturbation techniques in |HW99 1 have 
been used to obtain asymptotic expression for the price of vanilla options (call, put). 
Other cases have been derived using an asymptotic expansion of the heat kernel for 
short maturity (see |Lab05|). 

But Andersen and Andreasen in lAAOOII show that this sole assumption of lo- 
cal volatility is not compatible with empirical evidence (for instance, the post-crash 
of implied volatility of the S&P500 index). Hence, they derived a model with local 
volatility plus a jump process to fit the smile (we write it AA model). Their model 
may be seen as a perturbation of pure local volatility models. Of course, this is not 
the only alternative modelling The AA model fits some market data well (see cali- 
bration results in OAAOOl and those in this work), although we are aware that it does 
not work systematically nicely. In the following, we do not discuss the relevance of 
this model in specific situations. We simply focus on this model in order to illustrate 
our new approach for numerical pricing and fast calibration. For an analogous study 
on the time dependent Heston model, we refer to our work in preparation IIBGM08L 
Andersen and Andreasen lAAOOl calibrate their model by solving the equivalent for- 
ward PIDE. This sort of problem could be handled numerically using: an ADI-FFT 
scheme in [AAOO], a Finite Element Method in fMvPS041, an expUcit impUcit PIDE- 
FFT method for general Levy processes in |CV05 1 or Predictor Corrector methods to 
improve the accuracy of the PIDE in |BM06|. In the best case, all of these methods 
lead to a time of calibration of roughly one minute (see fAAOOl). Can we reduce this 
computational time? Is it possible to reach a time of calibration as short as the com- 
putational time of a closed formula such as Merton's IIMer76ll ? Our present research 
responds positively to the above questions. 

In order to handle even more general situations we consider, for the one dimen- 
sional underlying state process, the solution of the stochastic differential equation 
(SDE): 

dX,=^a{t,X,^)dWt + lJ.{t,X,^)dt + dJ,, Xq^xq. (I.I) 

For instance one may think of {X,), as the log asset price. Here {Wt)o<t<T is a standard 
real Brownian motion on a filtered probability space {Q,^, (^f)o<r<r,IP) with the 
usual assumption on the filtration {^t)o<t<T and {Jt)o<t<T is a compound Poisson 
process independent of (W,), defined by: 7, F, where Nt is a counting Pois- 

son process with constant jump intensity X and (y,),^^* are i.i.d. normal variables 
with mean rjj and volatility yj. Our main objective is to give an accurate analytic 
approximation of the expected payoff (or fair price of this option) 

miXr)) 

for a given terminal function h and for a fixed maturity T. 
The approximation can be applied to the following models: 



' for instance, see the book by Lewis | LewOO] on stochastic volatiUty models or the one by Gatheral 
IGat06l on models explaining the volatility surface. 
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Example 1.1 AA model on the log-asset. 

In this case, {Xt)t is the logarithm of the underlying asset, G{t,x) is its local volatility 

and ii{t,x) =X{\ — e''-'+T") — £_^jil jn order to guarantee the martingale property for 
[e^' )i . For a call exercised at maturity T, with strike K, h (jc) = -^o '-(u)du (^^jj {r{u)-q{u))du^x _ 
K)+ where r is the deterministic risk-free rate term and q is the deterministic divi- 
dend term. This model was derived in II AAOOII . In this work we mainly focus our 
discussion on this model. 

Example 1.2 Jump diffusion model on the asset. 

{Xt)t is the forward contract with maturity T, (7(f,x) is its volatility and pi{t,x) — 
—Xr\j. For a call exercised at maturity T, h{x) — e^ fo ''("^'^"(x — The primary 
focus of this model is the implied normal volatiUty instead of standard implied Black- 
Scholes volatility (Japanese markets in ||HKLW02i ) and it includes the presence of 
price jumps. 

Heuristics of our approximation and model proxy. In practice, at first glance, it is 
reasonable to think that (X,), (in the AA model) is approximated by a Merton model, 
where the coefficients jj. and CT only depend on time. We denote this proxy by (X/^); 
and it is defined by 

dXf ^a{t,XQ)dWt + jx{t,XQ)dt + dJt, X^^XQ. (1.2) 

This approximation can be justified by one of the following situations. 

i) The functions /i(-) and a{-) have small variations, which means that a(f,X,_) « 
a(f jjco) and analogously for /i. 

ii) The diffusion component is small (i.e. |/x|oo + |ct|oo is small) and the jump com- 
ponent as well (i.e. A(|77y| + jj) is small, meaning that the jump frequency or 
the jump size is small), which results in Xf w xq. This case is not equivalent to 
situation i) because the functions may be small and yet have large variations. 

Hi) Another obvious reason may be that the maturity T is small, leading to Xf k, xq. 

The heuristics i) and ii) are coherent with the parameter values taken in II AAOOl . 
When the three conditions are carried out at the same time, we expect our approx- 
imation to become even more accurate. Note also that no jump cases (A = 0) are 
allowed. The above qualitative features i) and ii) are encoded into quantitative con- 
stants Mo, M\ and Mj defined in ( 11.5b and will be discussed later in this work. The 
above heuristic rule implies that 

E(/z(X7')) = IE(/i(Xf )) + error. 

The term E(/i(X^)) is the price in the Merton proxy which is explicit (see Remark 
12.21 ). But this sole approximation is too rough to be sufficiently accurate. Our work 
consists of deriving correction terms for the above equality to attain a remarkably 
good approximation. 
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Smart expansion. To perform a rigorous analysis, we use a suitable parameteriza- 
tion w.rt. e e [0, 1]: 

dXf = e{a{t,Xl'_)dW, + jx{t,Xf-)dt + dJ,),X^ = xq, (1.3) 

so that X} = X,. We write 

g{e) = nh{^T)) (1-4) 

and our approach consists of expanding the price ( ll.4l l with respect to e. But the ac- 
curacy of the expansion is not related to e because the value of interest £ = 1 is not 
small. This is a significant difference as compared with singular perturbation tech- 
niques. Parameterization is just a tool to derive convenient representations. By using 
an asymptotic expansion in the context of small diffusions and small jumps (relative 
to the frequency or to the size), we can establish estimates of the derivatives. This 
allows us to make an explicit contribution at given order and to control the error. This 
is achieved by using the infinite dimensional analysis of Malliavin calculus. Here, 
we focus our analysis on the first termfl, for which we provide explicit formulas. 
We also give explicit upper bounds of the errors for general forms of ) and (t( ). 
However, the smaller the parameters and A(|77j| + jj) are, the smaller the 

maturity T is, or the smaller the derivatives of the functions o{-) and /i(-) w.r.t. the 
second variable are, the more accurate the expansion is. Given realistic parameters, 
the accuracy is indeed very good (less than 2bp in implied volatilities for various 
strikes and maturities). As a result of these expansions, we prove that the price (11.4b 
in our general model jl.lb equals the price in the Merton model plus a combination of 
Greeks (still in the Merton model). Hence, numerical evaluation of all these terms is 
straightforward, with a computational cost equivalent to the closed Merton formula. 
The residual terms (otherwise stated as error) are also estimated and their amplitudes 
depend on the smoothness of the payoff. We distinguish three cases: smooth, vanilla 
(call, put) and binary payoffs. In practice, the vanilla case is likely to be the most 
useful. 

This is our main contribution. Furthermore, from the approximation price we 
observe that one may obtain a volatility smile for short maturities (since we use the 
Merton model as a proxy) and a volatility skew for long maturities (due to local 
volatility function). 

Comparison with the literature. We refer in particular to Hagan et al. in I1HKLW021 

for the SABR model, to Fouque et al. in HFPSOOI for stochastic volatility models, or 
to Antonelli-Scarletti in I AS07I . In all these works, as opposed to our approach, a 
perturbation analysis w.r.t. the volatility, the mean reversion parameters, or the cor- 
relation, is performed and this leads to writing the price as a main term (essentially 
a Black-Scholes price) plus an integral of Greeks over maturities. In the time homo- 
geneous case, the authors successfully compute or approximate this integral, which 
strongly relies on PDE arguments. In our case, we do not approximate the underlying 
PDF (or the related operator) but owing to Malliavin calculus, we directly focus on 
the law of the random variable Xt given Xq = xq and not necessarily on the process 



^ in the former version of this work, terms at any order have been analyzed. 
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for any initial condition. Thus, we are able to handle time inhomogeneous coeffi- 
cients and jumps as well, without extra effort. This is a very significant difference 
from previous research. 



Outline of the paper. In the following, we present some notations and assumptions 
that will be used throughout the paper. Section|2]is aimed at presenting our methodol- 
ogy in an heuristic way to approximate the expected cost. Rigorous results are proved 
in Section|5] In Section[3] we derive financial modeling consequences from these for- 
mulas. These observations lead to justifying simplified choices of the local volatility 
(of the CEV type), to predict the form of all attainable smiles with their dynamics. In 
Section m we first give a methodology for implementing the approximation formula. 
Secondly, we show how to efficiently use our formula for calibrating the model using 
a relevant algorithm. Finally, we detail numerical applications in calibration for real 
market data using our simplified form of local volatility. In Section |5] we analyze the 
amplitude of the correction and error terms of the approximation formula; the anal- 
ysis depends on the kinds of payoff (smooth payoff in Theorem 15. II vanilla options 
in Theorem 15.21 binary options in Theorem l5.4b . In Appendix |6] we bring together 
useful results to make our "smart expansion" explicit. 



Notations used throughout the paper. 

Differentiation. If these derivatives have a meaning, we write: 

- yift '\x) — ^{t,x) for every function of two variables. 

- Xi^t = -^p-|£=o ■ These processes play a crucial role in the work that follows. 

- When there is no ambiguity, we simply write a, — a(f,xo),|ij = /z(f ,xo), cr/'' — 

The following definition is used to distinguish the payoff functions h. 

Definition 1.3 As per usual, we define '^q°°(R) as the space of real infinitely differ- 
entiable functions h with compact support (smooth payoffs). The sup-norm of the 
function h is denoted by \li\oo. We define as the space of functions with growth be- 
ing at most exponential. In other words, a function h belongs to Jif if \h{x) \ < cie'^^^kl 
for any x, for two constants ci and C2- 

The following notation provides a convenient representation of the correction terms. 

Definition 1.4 Greeks. Let Z be a random variable. Given a payoff function h, we 
define the Greek for the variable Z by the quantity (when it has a meaning) : 

Greek?(Z)^S^| ^. 

' V ^ g^, lx=0 

Given appropriate smoothness assumptions concerning h, one also has 

Greekf(Z) =E[/;(')(Z)]. 
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Assumptions. In order to get accurate approximations, we may assume that coeffi- 
cients (7 and jj, are smooth enough. 

- Assumption (R4). The functions (7(-) and /i(-) are continuously dijferentiable 
w.r.t. X up to order 4. In addition, these functions and their derivatives are uni- 
formly bounded. 

The functions and their derivatives could be piecewise continuous w.r.t. the time vari- 
able, without changing the following approximation formulas and the following error 
bounds. 

The assumption (/?4) seems to be restrictive because one requires a(-),/i(-) and 
their derivatives w.r.t. x to be bounded. On the one hand, this hypothesis is clearly 
too strong for us to use in the derivation of our smart expansion: indeed, the reader 
may check that polynomial growth conditions are sufficient for this purpose. On the 
other hand, assuming that the derivatives are bounded is much more convenient for 
explanation purposes. It enables us to state all our error estimates purely in terms of 
the following constants: 

Mi==maxi<K4(|cJ»|o„+|M''^|oe), 
Mo = maxo<K4(|cjW|oo+|Ai''^|oe), (1-5) 
Mj= \r]j\ + yj. 

Ml, Mq and Mj play complementary roles. 

a) The constant Mi is a measure of the norm of the derivatives (w.rt. x) of the objec- 
tive functions a{-) and /i(-). All our error estimates (see Theorems 15 . 1 H5 . 2115 .4b 
are linear w.r.t. Mi, which corroborates the proxy intuition explained in item i). 
The smaller the value of Mi is, the closer X and are, and as a result, approx- 
imation is increasingly accurate. At the limit Mi — 0, the initial model and the 
proxy coincide (Xt — Xf^) and our approximation formula becomes exact. 

b) The constants Mq and Mj also include estimates of the amplitudes of (7(-),/i(-) 
and of the jump components. All our error estimates also depend on powers of 
Mo and Mj. This mathematically justifies proxy intuition ii). The smaller Mq and 
Mj are, the better the resulting accuracy. 

In our next theorems, we also clarify the dependence of our estimates regarding jump 
frequency X and maturity T, because as these parameters decrease, the approximation 
becomes increasingly accurate. 

To perform the infinitesimal analysis, we rely on smoothness properties which are 
not provided by the payoff functions, but rather by the law of the underlying stochas- 
tic models (this is related to Malliavin calculus). The following ellipticity assumption 
on volatility combined with {R4) guarantees these smoothness properties. 

- Assumption (E). (7 does not vanish and for a positive constant Ce, one has 

1<^<C. 
where Gi„f = inf(,^,^)gi;+xK a{t,x). 
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We also need to separate our analysis according to payoff smoothness. We thus divide 
our analysis into three cases. 

- Assumption (Hi), h belongs to ^(J°(R). This case corresponds to smooth payoffs. 

- Assumption (H2). h is almost everywhere differentiable. In addition, h and h^^^ 
belong to M'. This case corresponds to vanilla options ( call, put). 

- Assumption (H^). h belongs to .y^f. This case includes binary options (digital). 



2 Smart Taylor Development 

In this section, we formally show how to replace the price E,{h{Xj)) by using that 
found in the Merton model E,{h{Xj)) with appropriate correction terms. Rigorous 
justification of the following expansions is postponed to Section|5] 

The initial trick of our smart expansion lies in the use of the parameterized process 
{Xf)t for e £ [0, 1], defined in il.3i . Under assumption (^4), almost surely for any t, 

is w.r.t e (see Theorem 2.3 in MFK85I ). If we put Xf^ = -^r^ we get 

dXl, ^o, {Xf_ )dW, + jX, {Xf_ )dt + dJ, 

+ eXl,_{al^\x^_ )dW, + lil^\xf_ )dt), Xf q = 0. 

From the definitions, Xi , = |£=o, c/'' = a^'^ {t,xo) and = ji^'^ (f ,xo), we easily 
get 

dXot=0, Xoo=xo, 

dXij=aidWt + fitdt + dJ,, Xi,o = 0, 

dX2j =2Xij_{ayUw,+f4^^dt), X2.0 = 0. 
Thus, the Merton model is obtained by the first order expansion of X^ at e = 0: 

Xoj+Xij —xq+Xij — Xj. 

We now use the Taylor formula twice: first, for Xj at the second order w.r.t e around 
e = 0, second for smooth function h at the first order w.r.t x around xq +X\ t — Xj . 
One gets: 

nK^'r)] =E[/;(xo + ^ + ...)] = ^[h{X^)]+^[h(') (X^)^] + .... 

Then, the price E[/i(X7-)] can be approximated by a summation of two terms : 

- '¥\h(Xj )]: The leading order which corresponds to the Merton price (BS price 
when A = 0) for the payoff h. 

- E[/i'''(X|^)^^]: The correction term which is made explicit in the next theorem. 
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Theorem 2.1 (Main approximation price formula). 

Suppose that the process data fulfills {R4) and (E) and that the pay off function fulfills 
one of the assumptions (Hi), {H2) or (H-i). Then 

3 3 

E[/i(Xr)] =E[/i(Xf )] + 2^ a;,rGreekf (Xf ) + ^ j3,- ^Greekf (Xf + Y') + Error, 

!=1 1=1 

(2.1) 

where 

rT rT 

aij = pL,{ y fii^^ds)dt, 

rT rT rT 

CCl.T = / {(yh / ^i^'-i^) +Mr( / C7,.(ji''t/i))c/f, 
Jo Jf Jt 

0:3. r = y asas^^ds)dt, 
Jo 

Pij =A / t{yji4^^ + r]ja,al^^)dt, 
Jo 

p3,T=^7j ta,a,^^^dt, 
Jo 

F' /s an independent copy of the variables 

/n addition, estimates for the error term Error /n the cases {Hi), {H2) and (i/3) are 
respectively given in Theorems \5.1\\5^ and \5.4\ 

To prove Theorem 12. II it remains to show that Wiji^^\Xj )^^^] is equal to the two 
summations of (12. 11 1. The reader familiar with Malliavin calculus for the computa- 
tions of Greeks (see ilFLL+99l . rGob04|, ...) may recognize in the expansion of 
E[/!(i)(X^)%^] the generic form of some derivatives (or Greeks) of E[/!(i)(Xf )], 
derivatives which are written as the expectation of h^^\Xj) multiplied by random 
weights. This is indeed our methodology to explicitly compute the correction terms 
in the formula ( 12.1b . 

Proof Define the new function G by G{x) = h{x + XQ + f^ jx,dt). One has: 



E&(i)(xf)]=E[^G(')(rMW,+7r)] 
2 2 Jo 



E[(/ Xi,,_((7/''t/W, + M/^'t/f))G"'( / a,dW,+JT)] 



Write (^[ ,)r for the continuous part of {X\ ,),. Using Lemma |6^ ( since Jj is indepen- 
dent of (Wi),g[0 7.]) andLeb{f e [QJ] : Xi., = = a.s. (see page 6 in ||Sat99l ), 
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one gets: 

E[^/z(i)(xf )] =E[( (\,ai'h{,dt)G^'-\r adW,+JT)] 
I Jo ' Jo 

+ E[(/ ,dt)G'^^\ o,dW,+JT)] 

Jo Jo 

+ E[( ^ a,Gl^h,dt)G^^\ r GidWr+Jr)] 
Jo Jo 

rT i-T 

+ E[(/ ii^^h,dt)G^^\ atdW,+JT)]. 
Jo Jo 

Apply Lemmas |6 . 3 1 and 16741 and use Definition 1 1.41 of Greeks to get the result. □ 

Remark 2.2 The above price approximation is a summation of three terms: 

1. W^hiXj )]: The leading order corresponding to the price when the functions (a,); 
and {lXt)t are deterministic. We know that in this case, there is a closed formula : 
either the Merton closed formula for call (put), or FFT tools for any other payoff 
because the characteristic function of Xj is explicit. For instance, the formula for 
a call in the Merton model (see LMer76J ) on the log asset is: 

where 

and BSCaU{S,K,T,v) is the Black-Scholes price for a call on an underlying St 
with initial condition So = S, volatility v, exercised at maturity T and strike K, 
where the risk-free rate and the dividend yield are set at 0%. 
2- LLi o:,,rGreek''(Z^): The volatility and drift correction term which depends on 
the first derivatives of jj. and a. This term can be computed as easily as the main 
term. 

3. ^f^i j3,, r Greek'' + Y'): The jump correction term which depends on the first 
derivatives of jx, a and on the jump parameters. Since Y' is also Gaussian and 
independent of Xj, the computation of these Greeks are similar to the previous 
ones, by adding to the mean Jq jj.tdt and variance a^dt the quantities rjj and 
7l 
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Remark 2.3 In the AA model on the log-asset, one has: 

2 Jo Jt Jo 

a2 r = - T / ctH <ys<yi^^ds)dt - A(e''-'+^ _ 1) / ta,al^^dt, 
2 Jo Jt Jo 

Pi.T = -Xr]j [ taay^dt, 
Jo 

Jo 

^3,7 =XYj / tatoi^^dt. 
Jo 

T f 1) 

Thus, the computation of these constants is simply reduced to that of /q tOtOt dt 
and a}{jj a,a^^^ds)dt. 

We note that we can perform higher order approximation formulas that remain ex- 
plicit. The only difference is that the number of random variables used as arguments 

for the Greeks will increase with each order, and it is within the set {Xi j +¥( H h 

Yl)i£n- We refer to I.Mir09i for higher order terms. 



3 Financial Modeling Consequences 

For simplicity, we consider the AA model on the log-asset (an analogous statement 
would be available for the jump diffusion model on the asset). 

The standard Gaussian framework as developed by Black-Scholes (1973) and 
Merton (1976) is realized by choosing a constant volatility function a{-) (the com- 
putation is still possible for a function dependent only on time). In order to arrive at 
a coherent, appropriate analysis and modelling for a fixed income market (without 
jump) Andersen and Andreasen IIAA02I take a parametric form for a: 

a(f,jc) = v(f)e('^W-')", (3.1) 

where v{t) the relative volatility function, j3 (f ) is a time-dependent constant elasticity 
of variance (CEV). PiterbargJ [PitOSJ uses the same form but applies it to Power 
Reverse Dual Currency swaps in order to handle the skew for the FX. 

Because of jU(f,x) — A(l — e''^+^) — ^-j^, the approximation formula ( 12. Il l 
depends only on a{t,xo),(y^-^\t,XQ),X,rij and yj. The volatility given in equation 
(13. Il l may generate all possible values of the following time-dependent functions 
(y{t,XQ) = v(f)e(/^W-')-^o and a'^^\t,xo) = (/3(f) - l)v{t)e^P'^'^-^^-'o, because it has 
two degrees of freedom v(f ) and /3 (f ). So this kind of volatility potentially creates all 

' If api, is the local volatility used in IPitOSI and L{t) = eio{i'(u)-q{ii))du ^ a{t,x) = apj,[t ,L,e^). 
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attainable prices in this class of models, and thus all attainable Black-Scholes smiles. 
This justifies interest in CEV-type volatility ( 13. 11 1. 

Attainable Black-Scholes smiles using the model. Can we predict the general 
form of the smiles generated by this model? 

- For short maturity: using our approach, the model is close to the Merton model 
related to Xj. Therefore, the shape of implied volatilities forms a smile centered 

on a point close to the money, which is on the left when Tjy + > (on the right 

when + < ). 

Formal Proof: Using the approximation formula, the correction terms are 0{T). 
So when T decreases to zero, the price converges to the Merton price. The second 
statement is easy to check. One can follow the approach of HGat02\\Mat00\l using 
characteristic functions, or can prove it directly using some derivations of the 
Merton formula liMer76]l . 

- For long maturity: the smile becomes a skew which is due to the local volatility 
function (because the smile for the Merton model flattens for long maturity). 

Smile Dynamics. The model has the Merton model as a good proxy. The implied 
volatilities for the Merton model are increasing and depend only on the ratio between 
the forward and the strike. Therefore, the smile should move in the same direction as 
the forward. 



4 Numerical Experiments 

In this section, we give details of the implementation for the approximation (12. Il l and 
illustrate the accuracy of our formula. After that, a generic bootstrap algorithm for 
calibration purposes is derived. Finally, a numerical application of this algorithm is 
applied to market data (currency options). 



4. 1 Numerical Implementation 

The case of time homogeneous parameters a;,(7/'',/Xr and /x/'^ gives us the coeffi- 
cients a and j3 exactly (see their expressions in Theorem l2.1l i. 
In addition, when these parameters are time-dependent, there are two cases. 

- Either the data are smooth. In this case, we use a Gauss-Legendre quadrature 
formula (see [PTVF92J ) for the calculation of the coefficients a and j3. 

- Or the data are piecewise constant. In this case, we can give explicit expressions 
of a and )3 in terms of the piecewise constant data. Let Tq = Q <Ti < ■ ■ ■ <T„ = T 
such that ai,Gj'^\lJ.t and /i/'' are constant at each interval ]7;,7/+i] and are equal 
respectively to aj;^, , Gjlj^,lJ.Ti+^ and /ij.''^ . Before giving the recursive formula, 
we need to introduce the following functions: coi^t = Jo o^ds, (02,t — Jo jJ-sds. 
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Proposition 4.1 Recursive formula. 

For piecewise constant coefficients, one has: 

, /-r ^\ (1) , (^'+1 ^ (1) 



a2,7;-+, =0!2^7;- + [T.+i - T^O^.^, «i,7;- + (^7;+, 4,+, cuj,?;) 
(7^-+i - 7i-)%^2 ,,(1) , ,, ^ ^(1) N 





=0:3.7; - 












^ii-^-^)„(i) 












=iS3.7; H 




















-7^')M7;+i. 



'1+1 7,-^| 1 



Proof According to Theorem l2.1l one has: 

«i,7;+i = IJ-ri J ds)dt + n,{ J n:, 'ds)dt 

=ai7;. + / /is 'ds)dt+ / / /Xi- 'ii5)t/f 

Jo 77; JTi Jt 

=ai.7;. + (/ llsds) jXtdt+ iJ.yds)dt 

JTi Jo JTi Jt 

=ai.T, + {Ti+i - T,)^^^l (02.Ti + ^^'^'~^'^\ t.^A'^, ■ 
The other terms are calculated analogously. □ 



4.2 Accuracy of the approximation 

Here, we give a short example of the performance of our method. The jump parame- 
ters have been set to: A = 30%, r]j = —8%, Yj = 35%. These parameters are not small, 
especially for the jump intensity A and the jump volatility Yj- The piecewise constant 
functions v and j3 defined in ( 13.1b are equal respectively at each interval of the form 
[2^, to 25% - / X 0.11% and 100% - / x 0.75%. The spot, the risk-free rate and 
the dividend yield are set respectively to 100,4% and 0%. 

We observe in the table below that the errors of implied Black-Scholes volatilities 
between our approximation and the price calculated using a PIDE method do not 
exceed 2 bp for a large range of strikes and maturities. The computational time of 
our formula is less than four milliseconds on a 2.6 GHz Pentium PC. The accuracy 
of our formula turns out to be excellent. 
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Table 4.1 EiTor in implied Black-Scholes volatilities (in bp=0.01%) between the approximation formula 
and the PIDE method expressed as a function of maturities in fractions of years and relative strikes. 



T/K 


70% 


85% 


100% 


120% 


150% 


3M 


0.02 


-0.03 


-0.92 


-0.07 


-0.12 


lY 


0.04 


0.06 


0.15 


-0.11 


0.01 


3Y 


0.22 


-0.23 


0.11 


0.41 


0.31 


5Y 


1.39 


1.06 


-0.01 


1.85 


1.76 



4.3 Calibration issues 

For this kind of model (AA model on the log asset or on the asset itself), calibration 
is still challenging as this model has no analytical formula. We can still perform a 
numerical calibration using the forward PIDE as explained in [AAOO|, but the time 
of calibration remains quite long (about one minute). With our approach, we can 
shorten the duration of calibration to less than one second, because our computation 
of the model price takes four milliseconds as previously mentioned. We achieve that 
by a simple bootstrapping algorithm using the path dependent formula. 

Bootstrap algorithm for piecewise data . Suppose that we want to fit option prices 
for n maturities Tq = < Ti < • • • < r„ and m strikes Ki, - ■■ , Km . First, we search 
the parameters X,r\j and Jj with best fit. At each interval ]7;_i, 7)], the data a, 0^^\ 

jj. and /x''' are constant, equal respectively to (Jj;, cjj.'', /Xj; and and depending 
on the vector Xi = {v{Ti),P{Ti)) (see formula [TTl i. Starting at / = 1, we express the 
coefficients ajj^ and Pj jj as a function of Xi, recursively using Proposition ( 14.1b . We 
apply a local minimization algorithm (for instance, the Levenberg-Marquardt as de- 
scribed in IPTVF921 ) in order to fit the implied volatilities for all strikes Ki, - ■ ■ ,Km 
at maturity 7; using our approximation ( 12.1b . Once the vector Xi is found, we go to 
the next step /+ 1, update a and j3 and compute Xi+i- 

This caUbration procedure is not completely safe. Sometimes we encounter insta- 
biUty problems. The final parameters depend on the initial guess. Moreover, there are 
many local minima. To avoid these problems, we could use a regularization method 
based on relative entropy (see [CT03 1), but these issues are not in direct relation with 
the accuracy of our formula. We think that the set of calibrated options (call/put) does 
not contain enough information on the future volatility to ensure a good calibration. 
Therefore, it is presumably worth including volatility options in the set of cahbrated 
instruments. This is a topic for further research. 

Calibration results. Here, we calibrate the EUR/USD exchange rate. The surface of 
implied Black-Scholes volatility is given in table |4.2| 

The jump parameters for the calibrated model are A = 1.21%, Tjy = —19.07% and 
Yj = 40.30%. The diffusion parameters v and j3 for the calibrated model are given 
in table 14.31 These values are realistic. The errors between the implied volatilities 
generated by the calibrated model and the market data are given in table 14.41 
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Table 4.2 Implied Black-Scholes volatilities for the EURAJSD rate expressed as a function of maturities 
in fractions of years and relative strikes. The spot is equal to 1.54. 



T/K 


92% 


96% 


100% 


108% 


6M 


10.82% 


10.65% 


10.53% 


10.56% 


lY 


10.84% 


10.70% 


10.63% 


10.66% 


1.5Y 


10.71% 


10.60% 


10.56% 


10.58% 


2Y 


10.60% 


10.48% 


10.46% 


10.47% 



Table 4.3 Calibrated values of the piecewise constant functions V and j3 . 



T 


V 


/3 


6M 


10.31% 


98.81% 


lY 


10.27% 


100% 


1.5Y 


9.90% 


100% 


2Y 


9.43% 


100% 



Table 4.4 Errors between implied Black-Scholes volatilities for the EUR/USD rate and those calculated 
within the calibrated model (in bp) expressed as a function of maturities in fractions of years and relative 
strikes. The spot is equal to 1.54. 



T/K 


92% 


96% 


100% 


108% 


6M 


-4 


3 


-1 


-3 


lY 


2 


1 





2 


1.5Y 


-1 


-3 


-2 


1 


2Y 


2 


-1 


1 


4 



The errors show that our model is a good model for the FX rate EUR/USD. Within 
our relevant algorithm, we are able to fit a 4 x 4 grid of quoted prices in less than 1 s. 

5 Error Analysis 

This section is devoted to the mathematical justification of Theorem 12. II and to the 
statement and proofs of upper bounds for the error term in (12. 11 1. For this, the analysis 
differs according to the payoff smoothness (smooth, vanilla or binary). We start with 
the smooth case (subsection l5.lt , which is less technical. Then, we handle the two 
other cases (call/put and binary options), which requires the use of Malliavin calculus. 

Throughout these computations, we aim at emphasizing the dependence of error 
upper bounds in terms of: the constants Mo, Mi and Mj defined in ( 11. 5t , the jump 
frequency X and the maturity T , in order to support the heuristic choice of the model 
proxy (see the discussion in the introduction). 

Additional notation. 
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- About floating constants and upper bounds. In the following statements and proofs, 
for the upper bounds we use numerous constants, that are not relabelled during 
the computations. We simply use the unique notation 

A <cB 

to assert that A < cB, where c is a positive constant depending on the model pa- 
rameters Mo, Ml , My, X,T,Ce (defined in assumption {E)) and on other universal 
constants. The constant c remains bounded when the model parameters go to 0, 
and it is uniform w.rt. the parameter e e [0, 1]. When informative, we make clear 
the dependence of upper bounds w.r.t. Mq, Mi, Mj, X and T. 

- Miscellaneous. As usual, the Lp-norm of a real random variable Z is denoted by 
||Z||p = [E|Z|''] In the proofs, the derivatives of the parameterized process 
are useful: they are defined by Xf^ = -g^- 

5.1 Error analysis for smooth payoff (under (Hi)) 

We begin our error analysis with the case of smooth payoff [h e '^/f^{M.)). 

Theorem 5.1 Error for smooth payoff. Assume that {R4) holds and that the payoff 
function h fulfills Assumption (Hi). Then the error term in Theorem \2.1\ satisfies the 
following estimate: 

|Error|<c sup \h^'\^{Mi^/f){{Mo^/ff+Mjy/xf). (5.1) 

;=1,2 

Let us briefly comment on the upper bound, making reference to the introduction. 
If the functions a(-) and /x(-) are only time dependent (Mi — 0), the approximation 
formula ( 12.11 ) is exact (the model and the proxy coincide). If they do not vary much 
w.rt. X {Ml is small), the accuracy is still good in view of ( 15.11) . If the coefficients 
(7 ( • ) , /X ( • ) and their derivatives and the jump size parameters are all small, the formula 
becomes very accurate. For instance, in a multiplicative case where (7(f ,x) — As{t,x), 
IJ.{t,x) — Am{t,x) and |Tjy| + 7y < 4 for a small parameter A, it readily follows that 
Mi,Mo,Mj== 0(4). Thus 

\Erroi\ = 0{A^T[Vf+VX]). 

Consequently, we may refer to the formula (12.1b in Theorem l2.1l as an approximation 
of order 2 w.rt. the amplitudes of the data (with error terms of order 3). 

These features arise similarly for the other examples of payoff smoothness. 

Proof It is divided into several steps. First, we write the SDEs satisfied by the three 
first derivatives of Xf^ w.r.t. e. Second, we give tight hp upper bounds on these deriva- 
tives. Finally, we combine these estimates with our smart expansion to complete the 
proof of Theorem |5.1| 
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Step 1. Differentiation ofX^. Under {R4), almost surely Xf is w.r.t e for any t 
(see Theorem 2.3 in IIFK85I ) and the derivatives are obtained by successive differen- 
tiations of the initial SDE il3i . Thus, direct computations lead to 

dXf, =a,{Xf_)dW, + n,{Xf_)dt + dJ, + eXf,_{al-^\xf_)dW, + ii^^\xf_)dt) , 

(5.2) 

dXl^[2Xf,_a^'\xf_) + e{Xf,_faPiXf_)]dW, 
+ [2Xl,_^^'\xf_ ) + e iXl,_ ) V/'' iXf_ )]dt 

+ eXl,_{a^^\xf_)dWr + ii^^\xf_)dt), (5.3) 
dXl,=[3Xl,_a,^'\xf_) + 3{Xf,_fa,^'\xf_) + 3eXt,_Xl_a^^^^ 

+ e{Xl,_faP\xf_)]dWr + [3Xl_^^'\xf_) + 3{Xf,_f^^^\xf_) 

+ 3eXl_Xl,_iiP {Xf_ ) + e{Xl,_ ) V/^'' {Xf_ )]dt 

+ eXl,_{a^^\x,'_)dWr + iij'^\x,'_)dt). (5.4) 

Their initial conditions are all equal to 0. Notice that unlike X^ and Xf, the processes 
X2 andXj are continuous. 

Step 2. Tight upper bounds. We aim at proving the following estimates for any p>2: 

E\Xl,\P <, {MoVTy +M'jlT, (5.5) 
mi,\P <, {MiVfy{{MoVfy +MPXT) , (5.6) 
mil" <c (MiVfriiMoVrfP+M^PXT), (5.7) 

uniformly for t <T. 

The existence of any moment is easy to establish, but here, we emphasize the de- 
pendence of the upper bounds w.r.t. the constants Mo,Mi,M/,A and T. Let us first 
prove the inequality (15.5b . From ( I5.2l i. apply Lemma |63] to the jump component and 
Burkholder-Davis-Gundy inequalities to the Brownian part, to deduce 

E|Xf J" <,f''/2-i ['E\a,{Xf)\Pds + tP-^ f'E\n,{Xf)\Pds+Mat 
Jo ' Jq ' 

+ tPl^-' ['E\Xf^a^^\xf)\Pds + tP-' ['E\X^,i^^''>{Xf)\Pds 
Jo ' ' Jo ' ■ 

<cTpI^mP+mPjXT + TpI^-^mP E\XfJPds. 

Using Gronwall's lemma, we easily complete the proof of (15.5b . For the second in- 
equality (15.6b . we proceed analogously and we obtain: 

E \P <cTP/^MP^ isupE\XfJP + sup E |Xf , | ^'^ ) . 

s<t s<t 
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Thus, plugging the estimate ( 15.51 ) into the previous inequality directly leads to (15.61 ). 
Now let us prove the inequality (15.7b . As before, apply BDG inequalities combined 
with Gronwall's lemma to obtain that 

IE|^3.,r <cr''/2Mf (supE|X|,|'^ + supE|Xf/'^ + supE\Xf, Xl,\P + supE|Xf,|^^'0. 

Use E|Xf^Z|J'' < i(E|Xf j2p + E|X|j2p) and the previous inequalities (|53E6l). 
Then bringing together different contributions easily leads to the required estimate 

(E3. 

Step 3. Completion of the proof. We follow the formal computations done at the 
beginning of Section|2] but more carefully. Let us introduce 

X2.T^ X^j{\-e)de, Xij= Xlj^—^de. (5.8) 
Jo ' Jo ■ i 

Then applications of Taylor expansions of at e = readily give these equaUties: 

Xt^Xj +X2J, Xt~Xj + -X2J +X^j 

where we have used Xj = xq +Xi j. Thus a second order Taylor expansion of h at 
point Xj writes 

E[h{X^)] ^E[hiX^ + ^+X,.T)] 

=E[h{X^)] +E[h^^\x^)^]+E[h^^\x^)X3,T] 

+ [\[h'^^\{l^v)X^ + vXT){X2Tf]{l-v)dv. 
Jo 

This proves that the Error term in (12. It for smooth payoff equals 

Error = E[h^^\x^)X^ j] + t E[h^^^{l - v)xf + vXt)Q12 r)^] (1 - v)dv. (5.9) 

Jo 

Then it readily follows that 

|Error| sup (E|X|^p)5 + sup E\X^t-\^. 

£G[0,1] ' £G[0.1] 

It is now straightforward to obtain Theorem 15. II by using estimates ( 15. 6115. 7b with 

p = 2. □ 

A careful inspection of the previous proof shows that assumption {Rt,) is sufficient to 
derive the error estimate ( 15. lb . 
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5.2 Error analysis for vanilla payoff (under (ffi)) 

This case has practical importance, because it includes call/put options. Regarding 
the error estimates related to Theorem 12. II we have paved the way with the case of 
smooth payoff. Nevertheless, there are some technical differences. The main one is 
that our previous proof represents the error in terms of the second derivative of the 
payoff, which is meaningless here. The additional ingredient is the Malliavin calculus 
integration by parts formula to avoid this second derivative appearing. We now state 
our main result when the payoff h is almost everywhere differentiable (with sub- 
exponential growth conditions). 

Theorem 5.2 Error for vanilla payoff. Assume that {R4) and [E) hold, and that 
the payoff function h fulfills Assumption {H2). Then the error term in Theorem 12.71 
satisfies to the following estimate: 

|Error| <, {p'\x^)\\2+ j\\h'^'\{l ~ v)Xf + vXT)hdv) 

x—{MiVf)({MoVff+MjVxf). (5.10) 

The shape of the upper bound regarding h is used for convenience in the proof. In 
view of the growth condition on /i''), the two first terms depending on /i''' are finite 
and uniformly bounded as Mq,Mi,Mj,X andr go to 0. 

Analogously to the smooth case (Theorem l5.1l l, the approximation error in ( 12.11) 
is of order 3 w.rt. the amplitudes of the model data, meaning that (12.1b is a second 
order approximation formula. 

Proof We split the proof into several steps. First, we assume that the payoff is smooth 
and we establish estimates that depend only on /i''^ the first derivative of h. For this, 
we need extra tools from Malliavin calculus, together with tight estimates on the 
Malliavin derivatives of the parameterized process. Then, we apply a density argu- 
ment to approximate h under {H2) by a sequence of smooth payoffs. 

Step 1. Malliavin calculus. For the usual Malliavin calculus on the Wiener space, we 
refer to Nualart IINua06l . But our case is slighty different because of jumps. However, 
in the following, our Malliavin differentiation is w.rt. the Brownian motion W and 
not w.r.t. the Poisson measure K. Hence formally, it is performed by leaving the jump 
component fixed, computing the Malliavin derivatives or integration by parts w.r.t. 
W, and then integrating out w.rt. the jumps. This principle has been formalized in 
several papers, for instance in IIBE08I Section 3. We briefly recall a few facts using 
their notations. 

The model jumps are associated with the Poisson measure K", with intensity gij^j^ {x)dx Xdt, 
where grij,yj is the Gaussian density on M with mean Tjj and variance The set of 
integer-valued measures on [0, r[ x R is denoted by Q^- For /(.) £ ^ — L2{[Q,T],M.), 
the Wiener stochastic integral /q l{t)dWt is denoted by Let denote the class 
of simple random variables of the form F = f{W{li), . . . ,W{In); k) where > 1, 
{Ii,...,In) G J^f'^, f : X 12^- R is bounded and infinitely differentiable w.r.t. its 
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first components (with bounded derivatives). We denote by D the Malliavin deriva- 
tive operator with respect to the Brownian motion. For F G .5^, it is defined as the 
^-valued random variable given by 

N 

D,F = ^ d,JiW{h), . . . ,WilNy, K)li(t). 

i=\ 

The operator D is closable as an operator from hp{£2) to Lp(f2,^), for any p>l- 
Its domain is denoted by D^ '' with respect to the norm || . ||i p given by = 

E\FY' \D,F\'^dt)Pl^. We can define the iteration of the operator D in such 

a way that for a smooth random variable F, the derivative D'^F is a random vari- 
able with values in As in the case k — \, the operator D'^ is closable from 
.y C Lp(i2) into Lp(i2; Jf®*), p>l. Its domain is denoted by D*^^ w.r.t. the norm 
\\P\\k.p = \&\F\P + Y.]=\^{\\D^F\\''^r,:i)YIP. With this construction, the operator D 
enjoys the same properties as the usual operator on the Wiener space (see |BE08l 
for more details). This justifies, in the case under study, the application of the usual 
results established without jumps (in particular the integration by parts formula and 
the related general Lp estimates, see the proof of Lemma |53] |. 

Step 2. Estimates of Malliavin derivatives. Under our regularity assumptions (^4), 
we know that for any f < T, any e G [0, 1] and any p > 1, we have G D'* '^ Xf^ G 
B^^P, G B^'P, X^^ G B^-P (see the arguments in IIBE08I ). Actually, we aim at 
proving the following tight estimates for any p>2: 



E\DrXf\P<,\a\P^, 




(5.11) 


EM^|'^<,|(T|£, 




(5.12) 






(5.13) 


E\Dl,^,X,r<c\cy\^M',P, 




(5.14) 


E\DrXl,\P<cMP, 




(5.15) 


E\Dl,Xlt\''<cM^MP, 




(5.16) 


E\DrXl,\P<cMP{{Mo^ 


^y'+M'j'XT), 


(5.17) 


E\dIxI,\p<,mPmP, 




(5.18) 


E\DrXl,\P<cMP{{Mo^ 




(5.19) 



uniformly in {r,s,t,u) G [0,7]'^ and e G [0,1]. Here again, the existence of any mo- 
ment is easy to establish and we will skip the details. We prefer to focus on the 
dependence of the upper bounds w.rt. Mq,M\ ,Mj,X and T. The bounds ( l5.16H5.18l - 
|5.19t are not used for vanilla payoffs, but only for binary ones. 
Proof of ^5J1\ . For r>t, DyX^ = 0. Now take r < f, in this case {DrXl^)r<t<T solves 
the following SDE (see IIBE08I ): 

/f 
D,X^_e{ai^\x^_)dW, + fi!,'\x^_)du), (5.20) 
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which defines a continuous process. Now, we proceed as in the proof of ( I5.5ll5.6ll5^ , 
combining BDG inequaUties and Gronwall's lemma. This gives 

^DrXf\P <, ¥.\Or{Xf_)\P + TpI^-^M'; j\\D,Xf,\Pdu <, \a\L 

and proves the announced inequahty. Besides, in light of (11.2b one has D,Xf^ — 
lr<t(yr, which directly gives ( I5.12l i. 

Proof of A5.13h Take for instance r < s <T, the other cases are handled in the same 
way. We have 

Dl,Xf ^eD,-X^_ai'\x^_) 

+ J'^ DX- £ {DrXf,_ {X^_ )dWu + Z),X„^_ Alp ) (x^ )du) 

which implies, in particular, that t ^ is continuous. It readily foUows that 

^Dljf\P <, Mf E|D,X/ 1" + TpI^MI sup V.\D,Xf,DrX^\P 

r<s<u<T 

<,MfE|D,X/|" + r"/2Mf sup (E\DsX^\^P + E\DrX^\^P)<c\a\'iMP 

r<s<u<T 

where we have used the Young inequality ab < ^ {a^ +b^) in the second line and 
(15.11b in the last inequality. The estimate (15.14b can be established in the same way. 
Proof of A5.15h We only consider r <t. Here one has 

+ J^D,XUd'\x^_) + eXt,,_a^^\x^_))dW., 
+ D,Xf,_ (Mi " {Xf,_ ) + eXL_ ' {XL ) )du 
+ D,-Xf „_ e (ai^^ {X^ )dW, + ^f/ ' {X^ )du) . 

It readily follows that 

nDrXl,\P <cnar{Xl) + eXl^_a^'\xL)\P + TPl^MP sup (E|D,-X„^|'' + E|D,X„%^, 

r<u<T 

Since a fixed time r is equal to a jump time with null probability and thanks to the 
Young inequality, we obtain 

nDrXl,\P <c \at+M'mx',^P + TPl'-M{ sup (E|D,-X,f |'' + E|D,X„^|2'^ + E|Xfj2''). 

r<u<T 

It remains to take advantage of the inequalities ( 15.5b and ( 15. lib , and to use |ct|oo < Mq 
and M\ < Mq to complete the proof of (15.15b . 

Proof of i5.16i5.17i5.1 8^5. 79l ). They can be proved similarly, with long and tedious 
computations. Since there is no extra difficulty, we will skip further details. 
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Step 3. Bounding the error using only h^^\ when h is smooth. We come back to the 
representation (I5.9l l for the error. The first term can be estimated using a Cauchy- 
Schwartz inequality and ( 15.71 ): 



\E[h^''\X^)X,j]\<,p'\x^)\\2 sup \\XIj\\2 

£G[0,1] 

<c\\h^'\x^)UM,VT){{MQVTf+M]^/Xf). 

This fits the required upper bound dS.lOt well, because Mq > Oinf. 

The second term in ( 15.91 ) requires a little extra work because of /z'^' . For this, we 
state a lemma, proof of which is given at the end. 

Lemma 5.3 Assume (E) and (R3). Let Z belong to np>iD^''\ For any v £ [0, 1], for 
k— 1,2, there exists a random variable Z],' in any Lp(p> I) such that for any function 

I G "^o^W' ^"^ 

E[/(*^) {vXt + (1 - v)Xf )Z] = E[1{vXt + (1 - v)X^)Zl]. 
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Moreover, one has \\Zl\\p <c 7S?Fi uniformly in v. 

Apply this Lemma with k = 1 and Z = (^2,r)^ defined in (I5.8I) . From the estimates 
(15.6115. 17l l. we readily obtain 

11^1 II 3 <. ^= . 

We have proved the upper bound ( 15.101 ). 



Step 4. Bounding the error under the sole assumption {Hj). So far, our error esti- 
mates depend on but they have been established for smooth payoffs h. It re- 
mains to justify that the error upper bound still holds for payoffs that are only almost 
everywhere differentiable (assumption (//2))- We argue by regularization, which is 
somewhat standard but a bit tricky here. We follow the proof of liGM05i . 

Denote by p the measure defined by lj^g{x)p{dx) = E(^(Xr)) + E(g(X^)) + 
E(g(Xf + Y')) + /o E{g{vXT + (1 - v)X^))dv. It is well known (see 0Rud66l for 
instance) that there exists a sequence {hf,)neN of smooth functions converging to h 
in L3(p) as well as its first derivative, as n goes to infinity. Thus, we can pass to 
the limit for E,{h„{XT)) and K{h„{Xj)). In view of (15.10b . we can also pass to the 
limit for the error bound. It remains to pass to the limit for the corrections terms, i.e. 
for the greeks Greek''" (Xj) and Greek|'" {Xj + Y'). To accomplish this, we represent 
them as E(/!„(X^)Z,) and E{h„{X^ + Y')Zi) using Lemma|53]with Z 1. Since Z; 
is in L3/2, we can pass to the limit as n ^ 00 to get E(/z(X^)Z,) — Greek'' (X^) and 

E(/i(Xf + y')Zi) =Greekf(Xf + y')- □ 



Proof of Lemma \5J\ Take = 1 or 2. 
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Step 1. Fy = vXt + (1 — v)Xy is a non degenerate random variable (in the Malliavin 
sense). Under [Ra), we know that Fy is in np>iD''^'''. One has to prove that 7/?, = 
Jq {DsFy)^ds is almost surely positive and its inverse is in any hp (p > 1). From the 
linear SDE ( 15.201 1 satisfied by {DsXt)s<t<T, we obtain 



Yf,. 



- [ + (1 - v)a^{xo)y-ds, 

Jo 

which clearly leads to our claim. Besides, for any p > 1, we derive 

\\rF,'\\p<c{cji„fVf)-\ 

Step 2. Integration by Parts formula. Using Proposition 2.1.4 and Proposition 1.5.6 
in IINua06l . one gets the existence of ZJ, in hp with 

<c ll7F„'ll^,2*p(2p+l)ll^^^'llt2*p(2p+l)ll^ll*:,p+^- 

Step 3: Upper bound of || 7^ ' far q > 2. On the one hand, using the 

inequahties (I5.11H5.12||5.13H5.14| |. we easily obtain 

ll^'^'vlk-., <. |cj|oo\/r. (5.21) 

On the other hand, with the same inequalities, we get sup,.<j.E|Dr7F,, <c T'^\a\^'M'^ 
and sup,.^,<7-E|D^,. 7/7,1'^ <c TP\(7\^'M^'' for any p > 2. Then, after some computa- 
tions, it follows that 

inf inf inf 

for any p >2. Finally using \a\oo < CEOinf (assumption (£)) combined with (15.211 ) 
and ( 15.221 1, we get 

This completes our proof. □ 



5.3 Error analysis for binary payoff (under (i/3)) 

For this kind of option, the payoff h is not necessarily smooth. We only assume that h 
is in J^. The results below are easy extensions of the vanilla options case, we leave 
the proof to the reader 

Theorem 5.4 Error for binary payoff. Assume that {R4) and (E) hold, and that the 
payoff function h fulfills Assumption (H^). Then the error term in Theorem \2.1\ satis- 
fies the following estimate: 

|Error| <, {\\h{X^)h + j\\h{{l ~ v)X^ + vXT)hdv) 

X (^ + ^)((Moyr)2+MiVAr). 
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Unlike the cases of smooth and vanilla payoff, for binary payoffs the approximation 
formula (12.1b is of first order w.rt. the amplitudes of the model data (with error terms 
of order 2). This is inherent to the lack of regularity of the payoff. 



6 Appendix 

6.1 Technical results related to explicit correction terms 

In this subsection, we bring together the results (and their proofs) which allow us to 
derive the explicit terms in the formula (12. 11 1. 

In the following, (ut) (resp. (v,) and (v,)) are square integrable and predictable (resp. 
deterministic) process and I is a smooth function with compact support. 



Lemma 6.1 For any continuous (or piecewise continuous) function f, any continu- 
ous semimartingale Z vanishing at t=0, one has: 




Proof This follows from the Ito formula applied to the product ( fsds)Zt . □ 
Lemma 6.2 One has: 

E[( r urdw,)i{ rvidw,)] = E[( r v,u,dt)i^^\ r vrdw,)]. 

Jo Jo Jo Jo 

In the case of deterministic u, it is equal to Jq VfUtdt Greekj [Jq VfdWt). 

Proof We first give the proof in a particular case when u and v are equal to 1 . By a 
usual integration by parts formula, one has: 

'E[1{Wt)Wt] = [ l{Vfx)Vfx^^dx= ( Tl'^^\VTx)^^dx^m[l^^\WT)]. 

For the general proof: apply the duality relationship of Malliavin calculus (see Lemma 
1.2.1 in |N ua06J ), identifying Ito's integral and Skorohod operator for adapted inte- 
grands. □ 

Lemma 6.3 Write {X'^^)t for the continuous part of {Xif)t. One has : 

E[( r v,X^tdt)l{ r v,dW,)] = r v,a,{ r v,ds)dtGreek\{ v,dWt) 
Jo ' Jo Jo Jt Jo 

+ Uti V,c/s)^^fGreeko( / VtdWt). 
Jo Jt Jo 
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Proof Applying first Lemma ISTI to /(f) = v, and Z, = one has: 
E[( f v,Xl,dt)l{ f v,dWt)]=n{ f { C Vsds)dXl,)l{ r VtdW,)] 

Jo ' Jo Jo Jr ' Jo 

= E[( r Vsds){a,dW, + ix,dt)l{ r vtdW,)] 

Jo Jr Jo 

= ( f VtOti r v,ds)dt)E[l'^^\ r v,dW,)] 
Jo Jr Jo 

+ {j\,{j\sds)dt)E[l{j\tdWt)], 
and we have used Lemma |6^ for the last equality. □ 



Lemma 6.4 One has: 

E[{ v,Jtdt)l{JT)]^X{nj r tVrdtGree\4){JT+Y') 
Jo Jo 

+ 7j / tv,dtGTeek[{JT + Y')), 
Jo 

such that Y' is an independent copy of the variables (y,),^^*. 
Proof Using the independence of increments for J, one has: 

EU r v,J,dt)l{JT)]= r v,E[J,l{JT-J,+J,)]dt^ r v,E[i{JT-J,)]dt. 

JO JO JO 

Using a conditioning ai-gument and since Yj is a Gaussian random variable, one 
has: 

k k 

l{x) = E[J,l{x + J,)] = ^ P(M = A:)E[ Yjl{x+Y. Yj)] 
kew j=\ i=\ 

k k 

= ^ F{Nr=k)k(rijE[l{x+Y^Yj)] + rjE[l^'\x+Y^Yj)]) 

kew j=i j=i 

k+l k+l 

= ^ ArP(M = k){iijE[l(x+ ^ Yj)] + 7,E[/(')(x+ ^ F,-)]) 
km 7=1 ;=i 

= Af (Tj7E[/ (x + /, + y')] + 7iE[/* ' H-^ + -^r + }'')]) , 

with Y' as in the lemma statement. □ 

6.2 Upper bound for compound Poisson process 

Lemma 6.5 The Lp norm (p> I) of the compound Poisson process at time t <T can 
be estimated as follows: 



E\J,\P <,MPXt. 
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Proof SetZj ~ (F/ — 77J )/);/. The random variables (Zy)y are i.i.d. Gaussian variables, 
with zero mean and unit variance. Then 

= 1 1 T?7 + 7jZj\ < \nj\Nt + yj\ £ Zj\ < Mj{Nt + I £ Zj\). 
7=1 ;=i 7=1 

Now it only remains to compute the p-th moment of A^, and Kt = lYfjLiZjl, which 
is considered a standard exercise. We give few details about the second term K,. 

First compute the characteristic function (p{u) =E(e"'^^=' ^) = exp(/lf(e"" — 1)). 
Then for an even integer p, one has Edjij Zj)P = ^{Kf) = iP(p<^P^ (0) = 0{^t). For 

odd values of p of the form p = 2k+ I, we apply the inequality ab < j(a^ + b'^) to 
write Kf < \ {Kj''+Kf''+^). The result then follows by using the estimates from the 
previous case {p even). □ 
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